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Abstract. In the present paper we describe the properties induced by disorder 
on an ultracold gas of Bosonic atoms loaded into a two-dimensional optical 
lattice with global confinement ensured by a parabolic potential. Our analysis 
is centered on the spatial distribution of the various phases, focusing particularly 
on the superfiuid properties of the system as a function of external parameters 
and disorder amplitude. In particular, it is shown how disorder can suppress 
superfluidity, while partially preserving the system coherence. 
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In the last years, the experimental results concerning confined ultracold atoms in 
optical lattices have attracted much theoretical interest from many different fields, 
ranging from condensed-matter physics to quantum information theory (see e.g. 
P, Q). The ability of tuning the fundamental physical parameters has provided 
an unrivaled tool to devise specific physical situations which, beyond their intrinsic 
physical interest, can be employed as important simulation tools [1]. In particular, 
the possibility to engineer a defect-free periodic potential appears as one of the 
most intriguing characteristics of these systems. Such "cleanness" has allowed the 
experimental realization of the Bose-Hubbard (BH) model and the observation of the 
supcrfluid-Mott quantum phase transition discussed by Fisher in his seminal paper 
where, in addition, the issue of the effect of disorder on the physical properties of 
the BH Hamiltonian was first addressed. Furthermore many experimental techniques 
such as laser speckle field [J] and the superposition of different optical lattices with 
incommensurate lattice constants d, S have recently substantiated the theoretical 
investigation of ultracold atoms in disordered optical lattices Q , revealing a rich 
scenario of newphysical situations such as the appearance of new phases (e.g. Bose- 
glass phase 0, 0j ) and superfluid (SF) percolation in d-dimensional (d > 1) lattices 

Mm. 

In this paper we will deal with the properties of a two-dimensional (2D) optical 
lattice where bosons are confined by an overall parabolic potential and subject to a 
random potential distribution. 

The experimental realization of 2D lattices is discussed in [13] , [3| and [3] • In 
particular the latter represents a direct experimental realization of the Hamiltonian 
discussed in the present paper, when disorder is absent. A valuable feature of the setup 
here considered is the possibility to investigate a wide range of physical situations and 
geometries according to different external parameter choices. While the parabolic 
potential confines bosons in a disk-like domain, the interplay of the other external 
parameters (hopping amplitude and chemical potential), allows one to realize, for 
example, Mott phases distributed in concentric shells, each with different filling, 
intercalated with SF shells. The ID counterpart of this scenario has been thoroughly 
studied in (lEj ]. 

Here we will focus particularly on the superfluid-phase spatial distribution in the 
presence of a random potential. We analyze first the situation where a (single) central 
disk-like Mott region is surrounded by a thin SF shell, and then the case when the 
system is fully SF. We show that, in the first case, the quasi-lD SF domain is strongly 
influenced by the presence of "impurities", leading to the drop of the SF fraction 
for a small increase of the disorder amplitude. In the second case, the genuine 2D 
geometry of the SF domain leads to the formation of circulating streams which do 
not disappear abruptly when the disorder amplitude is increased. In this respect our 
setup therefore appears to be an effective tool to investigate the combined effect of 
disorder and dimensionality on the SF properties of a bosonic system. 

The system considered is modeled by a BH Hamiltonian H = H] + Ht, with 

jj M M M 

Hi = ^ rii(ni - 1) - m rij, H T = -Ty]a\Aijaj + h.c. (1) 

i— 1 i—1 i,j 

where rii — a\a,i, a\ (<Zj) represents the creation (destruction) operator at site i, and 
the total boson number N = rij is a conserved quantity. The random on-site 
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effective chemical potential ^ is given by 

Hi = -n Q {x 2 t + y?) + Aj + (1 (2) 

accounting both for the noise distribution A (Aj = A-£j with £j uniformly distributed 
between and 1) and for the parabolic confinement £7o- The Hamiltonian parameters 
U and T represent the two-body interaction and the hopping amplitude between 
neighboring sites respectively, M is the total number of sites (in our simulations we 
have chosen M = 31 2 = 961) and the matrix Aij is the so-called adjacency matrix, 
which is zero if sites i and j are not connected, and one otherwise. 

Our analysis is based on the mean-field decoupling scheme already exploited 
by many authors for the analysis of the disordered BH model [H, [l(| [l6|. For an 



application of this scheme to a disordered system see 17|. Despite the well known 
shortcoming of not taking into account spatial correlations of quantum origin (see 
e.g. the discussion about spatial correlations in [ijj], and in particular Fig. 4 
therein), this mean- field decoupling scheme has proven the ability to disclose and 
qualitatively characterize the prominent physical features of Bosonic atoms loaded in 
optical lattices, expecially in dimension greater than one (for a detailed discussion on 
the MF scheme and its application see [19( and references therein) . Note however that 
we take into account the spatial dependence arising from the disordered+harmonic 
local potential, and hence — at least partially — the spatial correlations thereof. The 
essence of this approximation is the replacement of two-boson operators a|<Zj/ with 
single-site operators weighted by mean- field parameters an = (ip\ai\tp) (where a,; € C 
and \ip) represents the mean-field ground state). Namely 

ajfli' ~ a\aii + a*ai> — a*ai> with i S [1, . . . , N] (3) 

where o^s have to be determined by self-consistency standard procedure. 

Within this scheme, Hamiltonian H can thus be approximated by H ~ TLi + h 
where h = T^f^a^Aaae and parameters ji = ^2e=i Anai in the single-site 
Hamiltonian Tii = Un,i(ni — 1)/2 — in rii — T(7jaJ + 7*&i) take into account the original 
inter-site coupling. 

Most of the the fundamental physical properties of the system can be obtained 
from the values of calculated from Hamiltonian H alone where Tel. However, in 
order to properly describe the SF properties of the system it is well known that one 
has to go through same procedure with the introduction of the Peierls factors (20| . 
The Hamiltonian with the addition of the Peierls factors reads 

H = Hi + H T , H t = / '" /; ^.l,.,"., + h.c. . (4) 

The quantity OBij accounts for an infinitesimal velocity field v which can be 
considered as a "probe" to test the superfluidity properties of the system through the 
relation 8Bij = m/H f 3 v • dr 2l| where is the i-th site position. With the Peierls 

factors H becomes the Hamiltonian in the moving reference frame which is analogous 
to the Hamiltonian of a charged particle in presence of a magnetic field (see e.g. [2lJ). 
For the determination of the response of the system, it is necessary to take the limit 
9^0 and hence it is appropriate to consider the presence of the Peierls factors in 
perturbative terms with respect to the Hamiltonian H . Up to first order, it is possible 
to show that the mean-field parameters for the Hamiltonian 7i, corresponding through 
the mean-field site-decoupling approach to the Hamiltonian H, assume the following 
form ai = a" exp [—iOfa], where a° represent the (real) mean-field parameter related 
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to the Hamiltonian H., and the form of the phase term depends on the first order 
approximation in 6. 

In this paper we define the SF fraction (SFF) as the ratio between the current 
operator, as evaluated in Eq. ([7]) , and the maximum possible SF current between two 
sites, i.e. 

y..(ib\JiM) 

fs = = - 2T V^V^B^ (5) 

where we have employed the current operator defined as 

Jij = iT fexp[— iQBij]ajCLi — exp[i9Bij]ala,j\ (6) 

which describes the supercurrent in the rotating (lattice) frame. Recalling Eq. ([3]), 
the first-order approximation for the expectation value of the current operator can be 
written in the mean-field scheme as 

Jij = (V| J tj m * -2Ta°icPj6 [fa + Bij] . (7) 

So far, we have not explictly defined the structure of the site-dependent (through 
matrix Bij) Peierls factors. As suggested by the circular symmetry of the parabolic 
potential, we have considered by definition Bij an infinitesimal rotation around the 
paraboloid axis. The generated flow -observed in the rotating frame- must be 
construed as a response to this infinitesimal rotation. To interpret correctly Figs. 
[IJO and|H where the spatial distribution of Jij is depicted, it necessary to recall that 
supercurrents are depicted in the (rotating) lattice frame and hence, in the laboratory 
frame, the overall rigid-body rotation effect must be subtracted. The SFF at rest 
appears thus to be surrounded by a rotating rigid-body - in this case the non-SF 
fraction. 

In the ID homogeneous case f s vanishes iff = (the standard condition 
for Mott regime). In this case, the translational invariance of the system holds, 
entailing n, = rij = n, cq = a® = and <pi = <pj (overall phase coherence), 
thus Jij = iff a>o = 0. In 2D disordered case, conversely, Eq. shows that the 
condition f s = is ensured either by on — 0, like in the ID homogeneous case, or, 
more interestingly, by means of a phase rearrangement, which is determined by the 
constraint <j)j — <pi + B^ = 0, entailing the emergence of local clusters with aj =/= 0. 

It is worth noticing that, for a uniform (/i^ = fi) ID system, due to translational 
invariance, Eq. ([5]) reproduces the mean-field behavior of SFF as defined in the 
literature [II] 



The identification of / s id with a^/n deserves some comments. In 23] it was shown 
that the largest eigenvalue of the one-body density matrix pij — {ip\a\aj\il)) /N is a 
measure of the condensate fraction f c . Within the mean-field scheme, it is possible to 
prove that the ratio (af ) / {rii}- the space average (xi) (x = a, n) takes into account the 
lattice inhomogeneity due to the disorder-, can be related to the largest eigenvalue of 
the mean-field approximation of matrix = [((■0|n i |?/>) — |cti| 2 ) Sij + a*aj] \/N and 
thus provides an estimate of the condensate fraction f c . In this spirit, for a uniform 
system, Eq. © satisfies the condition f c ~ f s . From a more general point of view, 
since the condition f c ^ implies long range correlations, (af) / '(n,) can be viewed as 
a good coherence measure even for a disordered system. 
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In the numerical simulations we have analyzed two different situations. In the 
first case, we have evaluated the effect of noise on a configuration where a single central 
Mott domain, with rii = 1, is surrounded by a SF shell (Fig. [Tj), exhibiting a ring 
supercurrent. The parameters choice in this case has been performed in order to obtain 
a quasi-lD ring-shaped domain. In the laboratory frame the latter corresponds to inert 
matter decoupled from the lattice rotation owing to its SF character. This situation 
was explored in [3], in absence of disorder [24j. In comparison, as a second case, we 
consider a situation where a central disk-like SF region is present and, consequently, 
exhibits a 2D behavior. 




x/a 

Figure 1. Upper Panels. Expectation value for the local number operator 
(left) and value of (right). Lower Panel. Expectation value for the current 
operator Jij . Both configurations are considered in absence of disorder 



Case I. In the first case (T/U = 1.6 1(T 2 , ft /U = 6 1(T 3 and fx/U = 0.5 
corresponding to N to t = 199 total particles) the increase of the noise amplitude leads 
to the "fracture" of the ring with a t ^ and to the drop of the SFF, leading the 
rearrangement of the SF flow within the new potential pattern (Fig. [21). 

In the idealized situation of a perfect ID ring domain, the presence of an impurity 
cutting the ring is sufficient to determine a sudden drop to zero in the supercurrent 
induced by the velocity field. In the situation here depicted, due to the finite radial 
extent of the SF region such an impurity occupies an extended domain. In addition 
to this effect several vortices emerge which correspond to feeble supercurrents pinned 
around the sites with large (disorder-induced) local chemical potential fii (see Figf2J 
lower panel). It is worth noticing that in the laboratory frame, the extended impurity 
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Figure 2. Effect of increasing noise amplitude on the flux Jij (arb. units) 
distribution (left) and on the a; distribution (small panel, right). The figures 
considered represent noise amplitudes A/f2o = 60 (top), 67 (bottom). 



rotates together with the pinned vortices dragged by the lattice. 

The flow rearrangement in presence of disorder seems to be driven by the 
competition between the two terms in Eq. ([7]). In a situation where disorder is weak, 
the term Bij prevails, leading to a quasi-proportionality between the forcing term and 
the SF flow. On the other hand, in presence of large disorder, the phase pattern is 
adjusted so as to satisfy the condition <f)j — <\>i ~ —B^. In this case the parameter 
(af)/(n,-), a good SFF measure for a homogeneous ID system, while still being a 
measure of the local coherence of the system, can not be interpreted anymore as a 
measure of the SFF, since the presence of small clusters with on ^ will contribute to 
the overall (af)/(rii) value but not to the SFF. Notice that here by "small" we refer 
to a three-site cluster, in that, in a 2D square lattice 4 sites are needed to close a loop 
and thus to support a current. 

Fig. [3] shows that an increase of the disorder strength induces a sharp drop in 



Disordered 2D optical lattices 



7 



the SFF while leaving the quantum coherence of the system almost unchanged. Each 
data point is the result of the average over 10 realizations of the random potential. 
For sake of clarity, the errorbars are plotted every fifth data point. As we discussed 
above, the sharp drop in the SFF corresponds to the interruption of the quasi ID SF 
domain characterizing the system in the absence of disorder 

A qualitative argument for the determination of the disorder-amplitude range 
in which the ring starts to break is to consider that, in the absence of disorder, the 
ring is approximately 2-sites thick and with a radius of about 12 sites. Hence the 
potential energy difference between radially adjacent sites is about Qo ■ (l2 2 — ll 2 ). 
In order to have superfluidity, the hopping must be large enough to overcome the 
potential difference A c ~ Sl ■ (l2 2 — ll 2 ) , where A c can thus be viewed as the disorder 
amplitude corresponding to the drop in the SFF. As a consequence, if disorder-induced 
local potentials of up to the same amplitude are present, the SFF is not affected. 
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Figure 3. Comparison between the SFF (solid line) and the value of (ot?)/{ni) 
(dashed line) as a function of disorder amplitude. Error bars account for the 
dispersion of the numerical data around their mean value for 10 different disorder 
realizations. 



Case II. The second setup has been conceived in order to analyze the effects of 
disorder in a fully 2D SF domain. In this case, due to the more markedly SF character 
of the system, it has been necessary to consider much larger noise amplitudes to have 
observable effects on the SF distribution. In Fig. 2]we have sketched the values for Jij 
and a% in a situation exhibiting a strong SF character (T/U = 0.1, ilo/U — 3.6 10 -2 , 
fi/U = 0.6 with N to t — 313 ) for increasing values of the noise amplitude. Low values 
for Jij on the lattice center are direct consequence of the velocity pattern imposed on 
the lattice and must not be interpreted as the absence of superfluidity at the center 
of the trap. As previously discussed, in this case we have J7ij oc —Bij. 

For increasing values of the ratio between the noise and the parabola amplitude 
( A/fio ), as in Case I, the SFF, due to the rearrangement of the flow, has a slower 
decrease than in the previous situation, due to the much larger SFF - note the scale 
difference between Fig. [3] and Fig. [5] 

A prominent difference between Case I and Case II resides in the fact that the 
first one can be considered a quasi-lD case, and hence the "cutting" of the SF shell 
is analogous to the interruption of a ID chain. The second case represents, on the 
other hand, a genuine 2D situation and, as it is possible to see in Fig. [H the presence 
of noise leads to the formation of percolation patterns with nonzero circulation (Io| . 
The supercurrents in this case will follow closed and possibly merging paths encircling 
domains characterized by f s = but (af) / (rij) ^ 0. The latter shows then the survival 
of a local coherence character, connected to f c , even in the absence of supercurrents. 
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Figure 4. Effect of increasing noise amplitude on the flux Jij (arb. units) 
distribution (left) and on the a; distribution (small panel, right). The figures 
considered represent disorder amplitudes A/Qq = 200 (top), 600 (bottom). 



The condition f c ^ together with f s = 0, according to the scheme of Ref. [171 ] . 
supplies, in the thermodynamical limit, a possible characterization of a Bose-glass 
phase. 

In conclusion, we have presented a theoretical study of the behavior of Bosonic 
ultracold atoms in a 2D optical lattice with parabolic radial confinement with disorder. 
Our investigation has focused on the SF response both in the weakly and in the 
strongly SF regime. Within our approximation scheme, we have shown that, in the 2D 
disordered case, the absence of supercurrents is determined by two independent effects: 
the (usual) vanishing of the parameters o^, and the phase rearrangement implied by 
4>j — <f>i + Bij = 0. While the first effect can be related to the presence of a Mott phase 
in the thermodynamical limit, the second one, in the same limit, suggests the presence 
of a Bose-glass phase. In this framework, we have given a possible interpretation of 
the mechanism through which superfluidity is destroyed by disorder. Our numerical 
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Figure 5. Comparison between the SFF (continuous line) and the value of 
(a^)/(rii) (dashed line) as a function of noise amplitude. Error bars account 
for the dispersion of the numerical data around their mean value for 10 different 
disorder realizations. 



analysis evidences how, as a consequence of the phase rearrangement, the quasi-lD SF 
domain is destroyed by the appearance of "impurities" , with an ensuing drop of the 
SF fraction for a small increase of the disorder amplitude. Conversely, in Case II the 
genuine 2D nature of the SF domain leads to a multiple-stream SF flow determined by 
the increase of the disorder amplitude, with the consequent appearance of percolation 
patterns. In the laboratory frame, the strong-disorder sources in the rotating lattice 
have an overall site-dependent dragging effect on lattice bosons indicating the presence 
of non-superfluid regions locally surrounded by flows of nonzero circulation. 

The rich phenomenology here presented exhibits a complex interplay between 
parameters such as the hopping amplitude, the boson-boson interaction, the disorder 
amplitude and the parabola coefficient. In view of the intrinsic interest from the 
experimental point of view, this deserves a more systematic analysis we will perform 
elsewhere. 
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